Topology of the ground state of two interacting Bose-Einstein condensates 
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We investigate the spatial patterns of the ground state of two interacting Bose-Einstein con- 
densates. We consider the general case of two different atomic species (with different mass and 
in different hyperfine states) trapped in a magnetic potential whose eigenaxes can be tilted with 
respect to the vertical direction, giving rise to a non trivial gravitational sag. Despite the com- 
plicated geometry, we show that within the Thomas-Fermi approximations and upon appropriate 
coordinate transformations, the equations for the density distributions can be put in a very simple 
form. Starting from this expressions we give explicit rules to classify the different spatial topologies 
which can be produced, and we discuss how the behavior of the system is influenced by the inter- 
atomic scattering length. We also compare explicit examples with the full numeric Gross-Pitaevskii 
calculation. 

PACS numbers: 03.75.Fi, 05.30.Jp 



I. INTRODUCTION 

Bose-Einstein condensation of mixtures of different 
atomic species has recently been the subject of an in- 
tensive experimental and theoretical research [Q, |, |j| 

i 1 0, § |> 0, [|J- The first experimenta real- 
ization of a system of two interacting Bose-Einstein con- 
densates (BECs) has been obtained at JILA with a dou- 
ble condensate of 87 Rb in two different hyperfine states, 
\F, M F ) = 1 1 , — 1) and |2, 2) 0. This mixture was charac- 
terized by a partial overlap between the two condensates, 
in presence of a gravitational "sag" due to the different 
magnetic moment. Since then several other experiments 
have been performed with double condensates of rubid- 
ium |,|| and with spinor condensates of sodium in 
optical traps H. 

Motivated by these experiments and by the future pos- 
sibility of realizing other binary mixtures of interacting 
BECs, these systems have been extensively studied also 
from the theoretical point of view. Up to now only two 
particular cases have been addressed: (i) a system of two 
condensates with different mass in cylindrically symmet- 
ric potentials arranged concentrically || [?], P4 nil, and 
(ii) including a gravitational sag, but for condensates 
with the same mass (the JILA case) jl2| . 

In this paper we extend these studies by considering 
the very general case of two different atomic species, with 
different mass and in different hyperfine states, trapped 
in a magnetic potential whose eigenaxes can be tilted 
with respect to the direction of gravity. We show that, 
despite the complicated geometry, the ground-state of the 
system can be easily characterized within the Thomas- 
Fermi approximation valid for large numbers of atoms. 
We provide general formulas which allow to calculate the 



shape and the density distributions of the two BECs. Our 
results can be a useful tool to analyze future experiments. 

The paper is organized as follows. In Section II we dis- 
cuss the equations for the ground state of the system and 
show that upon an appropriate coordinate transforma- 
tion, they can be put in a very simple spherical form. In 
Section III we discuss the general features of the model, 
and we give an explicit algorithm to classify all the differ- 
ent topologies which can be constructed by varying the 
number of atoms and the inter-atomic scattering length. 
We also work out some example for the case 87 Rb and 
41 K, which is a promising system for the realization of a 
new binary mixture of BECs |l5|, [l6|]. Finally we com- 
pare the results against the numerical solution of the full 
Gross-Pitaevskii equations for the system, finding a good 
agreement. 



II. THE MODEL 

Let us consider a system of two Bose-Einstein conden- 
sates with mass rm and in the hyperfine state (Fi,Mpi), 
each containing Ni atoms (i = 1,2), confined in a mag- 
netic trap. The ground state of the system can be ob- 
tained by solving two coupled Gross-Pitaevskii equations 
for the condensate wave-functions ipi [|I3[ 
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with the normalization condition 
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The coupling constants Uij are given in terms of the scat- 
tering length a,ij by 1 12 
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where we used the fact that 0,12 
assume an, 022 > 0. 

The total potential experienced by each condensate is 
the sum of the gravitational potential and of a dipole 
magnetic potential U B (x) = fj,B9FiMpi\B(x)\ {gp% is the 
gyro-magnetic factor of the specie i) which we assume, 
as usual, to be harmonic 
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By defining K = (KiK 2 K a ) 1/3 , Xj = y/Kj/K, and 

Uoi = ^B9FiM Fl B 
u 2 = fJ-B9FiM Fl K/mi 

we can finally write Ub in the standard form 
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For what concerns the gravitational potential, here we 
consider the general case in which the vertical direction 
(the direction of gravity) is not aligned with any of the 
symmetry axis of the trap, but lies in the x — z symmetry 
plane, rotated by an angle 9. We include this possibility 
since in the experiments the trap confinement is generally 
weaker along the horizontal direction x, and therefore 
even a small angle can produce a large "horizontal sag" ; 
we will give explicit examples in the following section 
ifuj . The total potential is 

Ui(x) = U B l) (x) +m ig (x sin 9 + z cos 6) (11) 

By performing an appropriate transformation on the 
coordinates the potential U(x) can be put in a simpler 
form. These transformations amount to 

(i) a scaling by Xj (notice that the determinant of the 
transformation is equal to one) 



XjXj 



(12) 



in order to put {/b(x') in a spherically symmetric 
form; 

(ii) a rotation of an angle tp (x'j — > x") in order to align 
the z" axis with the vertical direction 
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The transformed potential reads (to simplify the nota- 
tions in the following we omit the apices, x't — > Xj) 



Ufa) = -m t Lof (r 2 + {z- z 0l ) 2 ) + U 0l (14) 
where we have defined r 2 — x 2 + y 2 , and 
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where the scaling factor / is given by 
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Then we perform a translation along z by zqi, defining 
dz = zq2 — zqi, and we express all quantities in dimen- 
sionless units, rescaling lengths by a^o = y/ti/ (muji) and 
energies by fiui. Finally, the expressions for the trapping 
potential that will be used in the rest of the paper are 



Vi(x) = EMx)-tfoi = \(r 2 + z 2 ) 



(18) 



V 2 (x) = C/ 2 (x) - U 02 = -V {r 2 + (z - dz) 2 ) (19) 



with 



dz 



m 2 L02 _ 9F2m F2 
miLof 9F\fnF\ 
Ig ( 1 1 



Ig ( m 2 



a ho \u% uj(J a ho u>{ \rjmi 



(20) 



- 1 . (21) 



To summarize, in this section we have shown that with 
suitable transformations the trapping potential for the 
two condensates can be reduced to a simple spherical 
form. This allows for a much easier investigation of the 
features of the interacting system, as will be discussed in 
the following section. 



III. THOMAS-FERMI APPROXIMATION 

For large number of atoms Ni the solution of Eqs. (|l|)- 
(||) can be derived in the so called Thomas-Fermi approx- 
imation which amounts to neglecting the kinetic terms 
V 2 ^- Therefore, by reabsorbing the values Uoi of the 
potentials on their minima in the definition of the chemi- 
cal potentials, ^i — Uoi —> Hi, the above equations become 

yi(x)+un|V'i| 2 + ui2|^2| 2 = Mi (22) 
V 2 {x) + u 2 i\^i\ 2 + U22H2? = H2 (23) 
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By defining 71 = 7721/7711, 72 
M11M22 — ^12: the solution of Eqs. 
lapping region take the form 
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Notice that in order to have overlap between ipi and 
-02 Eqs. (p7|)-(p8|) have both to be satisfied, that is 
both right members must be positive (l^i | 2 , \tp2\ 2 > 0). 
The overlapping region between the two condensates is 
therefore the intersection of the regions of space delim- 
ited by the spherical surfaces £j defined by the equation 
R 2 = r 2 + (z — z ci ) 2 , and identified by the sign of the 
coefficient on: for on > the region to be considered is 
the one inside the surface 5^, for a; < the one outside. 

In the regions where there is not overlap the wave func- 
tions take the usual form 

oil 2 = ^— (2/ii - r 2 - z 2 ) (35) 
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Analogously to the overlapping case, these solutions are 
defined in a region of space whose boundary is delimited 
by the surfaces Eoi of equation R\, L — r 2 + (z — z%) 2 , with 
2/ii, i?Q2 = 2/12/?/, z°i = an d z°2 — dz. 
Notice that in order to satisfy the continuity condition 
of the wave function ?/>j at the interface between the over- 
lapping and non-overlapping regions the wave function 
■001 must be connected to ip\ at the boundary defined by 
S 2 (where l^l 2 vanishes, but not l^il 2 ), and vice- versa. 



A. General considerations 

Even though a self consistent solution of the full prob- 
lem can be obtained only after having imposed the nor- 
malization of the wave functions, we can draw some gen- 
eral considerations by considering the role played by the 
determinant A and the coupling tj 12 . First of all we define 
the value of U12 at which the determinant A change sign, 
u = y/ui\U22- Then we notice that the behavior of the 
position of the centers along z, z C i, and the normalization 
factors Ui of the overlapping wave functions depends on 
two critical values 7i 12 = 777J11 and u\2 — "22/77, which 
define the poles of z C i and the zeros of a; (the latter have 
poles also for u\2 — It is not difficult to prove that 
one of this two values lies in the interval where A > 0, 
and the other outside. Therefore, to fix the hierarchy 
of the scattering lengths we choose the condensate 1 in 
order to satisfy the condition 77 2 > 7*11/7/22; with this 
choice the critical value lying in the interval of positive 
A is u* = rjUn. In Figs. |l| and ^| we show the position 
of the centers z C i and the normalization factors on as a 
function of W12. Notice that in correspondence of u* the 
center z C 2 of the surface £2 goes from —00 to +00, and 
the normalization factor ct2 becomes negative. Therefore 
for 77* < 77i2 < u the region where \ip2 \ 2 > is the one 
outside the surface E 2 (see also Figs. ||c,g). 

We can distinguish three cases: 

(i) 7ii2 < —7t, A < 0: no overlapping solution is al- 
lowed in this range. From Fig. || we see that 
both on are negative, and therefore it is not dif- 
ficult to prove that an overlapping region could be 
constructed only at the price of putting an hole in 
the condensate, where both ipi would be vanishing. 
This has obviously no physical meaning and in fact 
what actually happens is that when 7712 approaches 
—77+ the condensates eventually collapse |?], |l2| |. 

(ii) —77 < 7ti2 < tI, A > 0: in this range the two con- 
densate can coexist and overlap in some region of 
space if \dz\ < R w + R 2 o- We will discuss in detail 
the actual degree of overlap and its topology in the 
next section. 

(iii) 7ii2 > 7i, A < 0: in this case the strong mutual 
repulsion leads to a phase separation between the 
two condensates || ||. The actual shape of the 
interface is determined by the one at the critical 
value 77. Since for this value the overlap goes to zero 
(in the TF approximation) , if one further increases 
7ii2 the shape of the interface cannot change. 

Of course, if one retains the kinetic term in the GP 
equations this can in part affect the degree of overlap be- 
tween the two condensates. In particular the transition 
to the phase separation regime is not so sharp: the con- 
densates can have an appreciable overlap also for 7112 > u 
]To| , 0, [l^] . The effect of the kinetic energy is also to raise 
the critical value below which the system collapses. 
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FIG. 1: Plot of the rescaled position z c \jdz (continuous line) 
and Zci/dz (dashed line) of the centers of the "interacting" 
surfaces Ej as a function of the mutual coupling Uvi- 





1 «i 




-u 

"**\ ^ 

A<0 


1 u* ~~ s /' 

u 2i: 
1 

A>0 A<0 



FIG. 2: Plot of the normalization factors ot\ (continuous line) 
and C12 (dashed line) of the interacting wave functions as a 
function of 1112 (in arbitrary units). 



B. Topology of spatial configurations 

In this section we investigate the different configura- 
tions which can be obtained in the case (ii) discussed 
above (A > 0, — u < u\i < u). Before solving completely 
the system for some particular set of parameters, we give 
an overview of the different topologies that one can ob- 
tain. We again distinguish three cases, as shown in Fig. 



(I) "external overlap": This case can take place 
when the separation \dz\ between the centers is 
larger than the difference of the radii of the non- 
interactingprofiles, |JZoi — -R02 1 < \dz\ < R01+R02 
(see Figs. £3a-d). One can easily verify that all the 
four surfaces Ej and Eoi intersect on a circle per- 
pendicular to the plane in Fig. ||, passing for the 
points P and Q (shown as black dots). The overlap- 
ping region is the one contained between the sur- 
faces Ei and E2 (dashed lines in the figure) whose 




(d) S 



FIG. 3: Possible topologies for a binary mixture of two 
BECs. (1) "external overlap": U12 < (a), < U12 < u* 
(b), u* < U12 < u (c), and phase separation U12 = u (d); (2) 
"full overlap": U12 < (e), < U12 < u* (f); (3) "partial 
overlap": u* < U12 < u (g) and phase separation 1112 = u 
(h). Dark and light grey represent the regions occupied by 
the non-interacting condensates 1 and 2 respectively. The 
shaded area indicates the overlapping region. The bound- 
ary of these regions are delimited by the surfaces Eoi (non- 
interacting, continuous and dotted lines) and Ei (overlapping, 
dashed lines). 



actual shape depend on u\2, as shown in Fig. || for 
U12 < (a), < 1*12 < u* (b), u* < U12 < u (c), 
and M12 = u which is a case of phase separation 
(d). For smaller dz one obtains other configura- 
tions, which fall in the next two classes. 

(2) "full overlap": in this case, for \dz\ < \Rqi — -R02I, 
one of the two condensates is entirely contained into 
the other with whom it is fully overlapping. See 
Figs. ||e-f for U12 < and < 1*12 < u* respec- 
tively. Which of the two condensates lies in the 
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outer shell depends the actual value of the param- 
eters 

(3) "partial overlap": this is similar to the case (2), 
but now the overlap take place over a shell which 
separates the inner core containing the condensate 
1, and the outer region with the condensate 2 (Fig. 
||g). Notice that according to the above discussion 
this configuration is possible only for u* < u±2 < u 
where the sign of 012 is negative. This is a nec- 
essary but not sufficient condition, since also the 
condition R2 < R\ must be satisfied. In this range 
of U12 (where 012 < 0) another possible solution is 
i?| < 0, which leads to the case (2). By further in- 
creasing U12 to the critical value u one again obtain 
a configuration of phase separation (Fig. ||h). 

Having determined the possible configurations of the 
system, we are now ready to solve any particular problem 
by imposing the normalization condition (||). To do this 
one has to write the normalization integrals for each of 
the possible profiles in Fig. ||, and then solve Eq. (||) in 
order to find the chemical potentials /i.; as a function of 
the atom numbers iVj. The analytic expressions for these 
integrals are given in Appendix |A|. These are polynomial 
functions of fractional powers in the chemical potentials 
fii , and in general Eq. (|3|) does not admit analytical solu- 
tions. Therefore the relation between \ii and Ni must be 
inverted numerically (which is nevertheless a much easier 
task than solving numerically the full Gross-Pitaevskii 
problem). For the special case of phase separation the 
two normalization equations can be decoupled (by using 
the fact that R\ — R2 and z c i = z C 2 for U12 = u), and 
solved analytically for dz = 0. 

Notice that in general (except for some particular case, 
e.g. dz = 0) it is not possible to know a priori which of 
the various configurations in Fig. |^ applies: one has to 
solve Eq. (|J) for all the possible configurations, and then 
choose the one which gives a self-consistent solution. 

In summary the ground state configuration for a par- 
ticular system can be found in three steps: 

i) choose the normalization integrals which applies 
to the possible profiles in Fig. || for a given U12, 
and determine /ij(iVj) by solving Eq. (||) self- 
consistently; 

ii) identify the overlapping region by plotting the "in- 
teracting" surfaces £j; 

iii) determine the non-interacting region for each con- 
densate by using the "non-interacting" surfaces 
Eoi, and the continuity of the wave functions (see 
Fig. |). 

We also remind that for the very special case dz = 
there are also symmetry breaking solutions, not included 
in the present analysis, which could be energetically fa- 
vorable M. 



C. Examples 

To give some explicit example we now consider two 
condensates of 87 Rb and 41 K, which is a promising sys- 
tem for the realization of a new binary mixture of BECs 
[PI jlj). We will classify some possible configurations 
which can be obtained by varying the number of atoms 
in each condensate, for different values of the inter-atomic 
scattering length, which is considered here as a tunable 
parameter The results, valid in the Thomas-Fermi 
(TF) approximation, will be compared with the numeric 
solution of the full-3D Gross-Pitaevskii equations (GPE), 
found using a steepest descent method [ |13[ . 

We start by considering a case of "external overlap" 
with both condensates in the hyperfine level |2,2) (77 — 
1). The scattering lengths are — 99 ciq and ax = 
60 ao, a being the Bohr radius |0. As trap frequencies 

j™ = 250 Hz, with an 



we use uj 



Rb _ 



16 Hz, w, 



Rb 
V 



angle of rotation 9 = 0.035 (we retain these values for all 
the cases analyzed in this section). With this choice the 
reduced coupling constant 



u Rb ,Rb = 0.0611 , u KiK = 0.0785 



(37) 



and therefore, according to the above discussion, we iden- 
tify the condensates 1 with 87 Rb, and the condensates 2 
with 41 K. We choose a case attractive interaction between 
the two condensates, a\2 = —55 ao (1*12 = —0.0530), with 
N Rb = 5 • 10 4 and N K = 2 • 10 4 . To visualize the role 
of the scaling and rotation transformations, in Fig. |] we 
show the TF profiles of the two condensates in the x — z 
plane, in rescaled (left) and natural coordinates (the co- 
ordinate axes correspond to the trap eigenaxes; right). 
The profiles in natural units can be easily obtained by 
performing the inverse transformation of those in Eqs. 
@ and Q. We will use this system of coordinates to 
show all the following figures. Notice that despite the 
small rotation angle 9, (the direction of gravity, repre- 
sented by a dotted line in the right picture of Fig. ||, is 
almost indistinguishable from the z axis on the scale of 
the figure) the misalignment in the direction of gravity 
produces a relatively large horizontal sag in the x direc- 
tion where the trap confinement is weak. 
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FIG. 4: TF profiles of the two condensates (non-interacting: 
continuous; overlapping: dashed) in rescaled (left) and natu- 
ral (right, in units of aho) coordinates [x horizontal, z vertical) 
for a case of "external overlap": Nm, = 5 ■ 10 4 , Nk = 2 • 10 4 , 
Oi2 = —55 ao. 

In Fig. [5] we compare the TF profiles with the contour 
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plot of the two densities, as found from the full GPE 
solution. For clarity each condensate is plotted sepa- 
rately, and compared with the TF profiles which define 
the boundary of the non-interacting (continuous lines) or 
overlapping phases (dashed lines), as defined in Fig. [|. 
The outer contour line for each condensate correspond to 
10% of its peak density (for y = 0). 
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FIG. 6: Density contours of the GPE solutions in the x — 
z plane for the 87 Rb (bottom) and 41 K (top) condensates, 
for a case of "partial overlap": Njo, — 2 ■ 10 4 , Nk = 2 ■ 
10 J , ai2 = 67 ao. Each condensate is compared with the TF 
profiles which define the boundary of the non-interacting or 
overlapping phases, as defined in Fig. ^g. Lengths are given 
in units of a^o- 



FIG. 5: Density contours of the GPE solutions in the x — z 
plane for the 87 Rb (bottom) and 41 K (top) condensates. Each 
condensate is compared with the TF profiles which define 
the boundary of the non-interacting or overlapping phases, 
as defined in Fig. ^a. This a case of "external overlap" with 
attractive interaction between the two condensates, ara = 
—55 ao, and Nm> = 5 ■ 10 4 , Nk = 2 ■ 10 4 . Lengths are given 
in units of a;, . 

Then we consider two examples for a system where the 
87 Rb condensate is in the hyperfine level |2,2) and the 
41 K condensate in |2, 1) (we use again clk.k = 60 ao). In 
this case 77 = 0.5. In Fig. || we show a case of "partial 
overlap" , obtained by fixing the inter-atomic scattering 
length to a±2 = 67 ao (W12 = 0.0645) and the number of 
atoms to N Rb = 2 • 10 4 and N K = 2 • 10 5 . Notice that 
when both condensates are in the |2, 2) level the spatial 
separation between the two is too large to allow for a 
configuration of "partial overlap" for reasonable values 
of the trap frequencies (in principle one could reduce the 
separation by strongly increasing the confinement in the 
direction of gravity). 

Finally, in Fig. ffl we show a case of "full overlap" for 
N Rb = 5 • 10 5 , N K = 1 • 10 4 , and a u = 20 a , giving 
uu = 0.0193. 

From the examples considered here we see that, al- 
though the full solution of the GPE is required for a 
precise determination of the actual degree of overlap be- 
tween the two condensates, the TF approximation well 



captures the basic topology of the ground state configu- 
rations. Therefor, due to its simplicity, the TF method 
presented here can be a useful tool to characterize the 
ground state structure of a binary mixture of BECs also 
in presence of a non trivial geometry. 

We conclude this section by noting that we have also 
verified that our method well reproduces the results al- 
ready studied in literature in case of simpler geometries 



IV. CONCLUSIONS 

We have presented a general method to classify the 
ground-state of a binary mixture of Bose-Einstein con- 
densates. We have considered the general case of of two 
different atomic species, with different mass and in dif- 
ferent hyperfine states, trapped in a magnetic potential. 
We have explicitly included the possibility of a non triv- 
ial gravitational sag, when the direction of gravity is not 
aligned with any of the trap eigenaxes, since even a small 
misalignment can produce a large "horizontal sag" . We 
have shown that, within the Thomas- Fermi approxima- 
tions and by performing a suitable coordinate transfor- 
mation, the equations for the density distributions can 
be put in a simple spherical form. We have given explicit 
rules to classify the different spatial topologies which can 
be produced, and we have discussed how the behavior of 
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FIG. 7: Density contours of the GPE solutions in the x — z 
plane for the 87 Rb (bottom) and 41 K (top) condensates, for 
a case of "full overlap": N Rb = 5 • 10 5 , N K = 1 • 10 4 , 
di2 = 20 ao- Each condensate is compared with the TF 
profiles which define the boundary of the non-interacting or 
overlapping phases, as defined in Fig. ^f. Lengths are given 
in units of a;, D . 



came aware of a very recent preprint related to this sub- 
ject |§. 



APPENDIX A: NORMALIZATION INTEGRALS 

In this appendix we give the explicit expressions for 
the integrals which enter the normalization condition <^). 
We distinguish two general cases: (i) "internal overlap" 
one of the two condensates is entirely contained into the 
other (Fig. |3|, right column), and (ii) "external overlap", 
(Fig. ||, left column). 

In both cases the normalization condition can be writ- 
ten as a sum of integrals of a generic density 

= a {Rl ~ r 2 - (z - z c f) (Al) 

over an appropriate portions of spherical domain whose 
boundary is given by the surface 

R\ = r 2 + (z - z A f . (A2) 

In the following subsections we consider explicitly the two 
cases. 



1. Internal overlap 



the system is influenced by the inter-atomic interaction. 

We have also provided explicit examples, and com- 
pared the results with the full numeric Gross-Pitaevskii 
calculation, finding a good agreement. 

The results presented in this paper might be useful for 
analyzing future experiments where new combinations of 
binary condensates are likely to be produced |ll| [ll| . 



In this case the normalization condition can be im- 
posed by using a combination of integrals over spherical 
domains. The generic form is 



I T (a, z c , R c , za, Ra) 



inaRi 



-R„ 



1 



(z c - z A f - \r 



(A3) 
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From this expression one also recover the value of the 
integral for the non-interacting case 



&7T 

In(Hi) = -T7 aR l 
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(A4) 



Note added After having completed this work we be- 



(I n (jii) = 4 7 r(2/i l ) 5/2 /u l4 for V = 1). 

B y using appropriate combinations of the integral 
( [A3] ) ,_the normalization condition for the case shown in 
Fig. |g reads 



J 



iVi 



JVo 



Ii{a!, Zci, R%(ij,i, fj, 2 ), z c i, Riini, - Ii(ai, z c i, Ri(nx, z C 2, R2{^i, 1^2)) 
+I I (a 01 ,0, splp^fq, z c2 , R 2 {^i,^2)) 
I n (n 2 ) - Ii(a 02 ,dz, v^^/f?, z c i, Ri{nx, 112) 



(A5) 



+Ii(a 2 , z c2 , R 2 (m, fx 2 ), z c i, R%(ni, fi 2 )) - Ii(a2,z c2 ,R 2 (ni,{i 2 ),z c2 ,R 2 (ni,iJ, 2 )) 



(A6) 



where we have indicated the explicit dependence on /ii and /i 2 - The cases in Figs. |^,f,h can be constructed in a 
similar way. 



2. External overlap 



These are the configurations shown in Fig. ^p-d. In this case the master integral can be written as the integral over 
a convex domain delimited by two spherical surfaces (A and B the upper and lower ones along the z axis respectively) 

I E (a, z c , R c , z A , R A , z b , Rb) = tto [R 2 a (R 2 c - 0.5R 2 A - (z c - z A ) 2 ) (R A ~ z(A, B)) 

-R 2 A (z c - za){R\ - z 2 {A,B)) - ± (R 2 C - (z c - z A ) 2 ) (R 3 a ~ z 3 (A,B)) 



+ \(z c 



z A ){R\ - z 4 (A,B)) + -(R A - z 5 (A,B)) 



+ {za *-> zb,Ra 



-Rb) 



with 



z(A,B) = f^ 



g| + (fg - za) 2 
2(z A - z B ) 



(A7) 



(A8) 



By assuming a configuration where the condensate 1 has a lower position along z (as in Fig. |^), the normalization 
condition for the cases with m 2 < u* shown in Fig. ^a,b is 



Nl = In(fJ-l) - lE(aoi,0, \/%Ll,0, V^Ul, Z C 2,R2{H1,^2)) 

+I E (ai, z cl , Ri(ni, /i 2 ), z c i, Ri(ni, Ma), z c2 , R 2 (pi, Ma)) 
N 2 = ^n(Ma) - Ie{<*Q2, dz, y/2/j, 2 /r), Zd, Ri(pu M2)) dz, V^Mi) 
+Ie(oc2, z C 2, R 2 (hi, M2), z c i, Ri(Hi, Ma), z c2 , R 2 (pi, Ma)) 

I 



(A9) 
(A10) 



In an analogous way one can construct the appropriate ( "external overlap" ) , by considering the appropriate corn- 
normalization condition for all other cases in this class bination of integral of the form (A7) over convex domains. 
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